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Abstract: We consider various sufficiently nonlinear sigma mod- 
els for nematic liquid crystal ordering of RP 1 ^^ 1 type and of lattice 
gauge type with continuous symmetries. We rigorously show that 
they exhibit a first-order transition in the temperature. The re- 
sult holds in dimension 2 or more for the RP 1 ^^ 1 models and in 
dimension 3 or more for the lattice gauge models. In the two- 
dimensional case our results clarify and solve a recent controversy 
about the possibility of such transitions. For lattice gauge models 
our methods provide the first proof of a first-order transition in a 
model with a continuous gauge symmetry. 



In jTU], it was shown that a class of ferromagnetic N- vector models with a 
sufficiently nonlinear nearest-neighbour interaction - meaning that the near- 
est neighbour interaction has the shape of a deep and narrow well - show a 
first-order transition in temperature. An example of a model which could be 
treated is given by the Hamiltonian 



with p large. These results confirmed earlier numerical work jHlEl- The main 
ingredient of the proof was a similarity between such models and high-g Potts 
models, which allows one to adapt proofs for Potts models, first developed 
in ^7j, and based on |SJ, to prove first-order transitions in the temperature 
parameter for these iV-vector models. 

In this paper we extend our analysis to include RP 1 ^^ 1 (liquid-crystal) 
models (such as were first introduced by Lasher and Lebwohl [2HJ l22j) and 
lattice gauge models. We can then employ techniques for showing the exis- 
tence of first-order transitions, such as have been used for Potts ferromagnets 
in d at least 2, as well as for Potts lattice gauge models in d at least 3 [T7] . 
to conclude that the corresponding nonlinear liquid-crystal and lattice gauge 
models, (with either abelian or non-abelian symmetries) have 1st order tran- 
sitions. 

The standard iV-vector models are either believed or sometimes rigor- 
ously known to have 2nd order transitions in d = 3 or higher, a " Kosterlitz- 
Thouless" transition in d = 2, iV = 2, and no transition for d = 2 and higher 
N pp. In the XY-model (N = 2) for either d = 2 or high d these results are 
rigorous, for the other models there is a consensus based on both numerics 
and heuristic arguments. 

In contrast, for the standard versions of the liquid crystal and lattice 
gauge models, as well as for very non-linear cr-models, both numerics and high 
temperature series suggested the existence of 1st order transitions, despite 
some theoretical analyses originally either suggesting 2nd order transitions, 
no transitions at all, or Kosterlitz-Thouless type transitions. For some of 
this literature, see e.g. [2312312211111121122112312111211 and references 
therein. Moreover, in the limit where N approaches infinity (the spherical 
limit) 1st order transitions were found, in dimension 2 or more f§l EHl EH] ■ 
This spherical limit result also holds for our nonlinear interactions in the 
ferromagnetic case jjj- Whether such a first-order transition can also occur 
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for finite iV larger than 3 in d = 2, or whether it might be an artefact of 
the spherical limit has for a long time been a matter of controversy (see for 
example [HU EH])- In fact, Sokal and Starinets described the existence of 
such a first-order transition as a "pathology". Our result finally settles this 
question: first-order phase transitions for models with continuous symmetry 
in d = 2 can occur, despite the conjecture to the contrary of |3T]- Our results 
are essentially in agreement with the analysis of in contrast to what 
was suggested in most earlier analyses, the symmetry or the low-temperature 
properties of the model do not play a role of any great importance, and in 
fact for our nonlinear choice of interaction the spin-dimensionality N does 
not need to be large and can be as small as 2. The lattice dependence of the 
phenomenon found in JS2J seems somewhat of an artifact, however, which 
disappears if one varies the nonlinearity parameter. 

The fact that our proofs are insensitive to the nature of the phases be- 
tween which the transition takes place implies that, similarly to the ferro- 
magnetic cases, one might have in the liquid- crystal models a transition be- 
tween a disordered high-temperature phase and either a nematically-ordered, 
a Kosterlitz-Touless or a disordered phase at low temperatures. Similarly 
one might find a transition either between a confining and a nonconfining - 
Coulomb- like - phase or between two confining phases in the lattice gauge 
models. Which one occurs in a particular case should depend on dimension 
and/or symmetry of the system. 

In particular, we emphasize that our proof is also insensitive as to whether 
the symmetry group of the lattice gauge model is abelian - in which case 
it is expected that in 4 dimensions a transition between a confined and a 
Coulomb-like phase occurs ^HJ ^] -, or nonabelian, in which case both 
states are expected to confining, (this is also expected in general in d — 3). 
For a heterodox discussion on the difference between what is to be expected 
in abelian and nonabelian models, including some history of this problem cf 

We consider a lattice Z d , and either spin models in which the random 
variables live on the sites, or lattice gauge models where the variables live 
on the bonds (or links) of the lattice. The parameters of our models are the 
spatial dimension d, the spin-dimension N, and the nonlinearity parameter 
p. The "standard" versions of the models are obtained by taking p — 1. 

For ferromagnetic models the variables are iV-component unit vectors 
living on the iV-sphere. We will discuss the argument in the 2-component 
case, in the general case the proof is essentially the same. The ferromagnetic 
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models treated in were given by 



H = -J Y, 1 + COS f ~ <P]) ) • (1) 

For liquid crystal RP 1 ^ -1 models either one considers variables - usually 
denoted n« - which live on the projective manifold obtained by identifying a 
point on the iV-sphere with its reflection through the origin, or equivalently 
one can consider ordinary spins on the iV-sphere, and divide out this "local 
gauge symmetry" afterwards. The last approach is the route we will pursue, 
as it allows us to literally transcribe the proof of [TOj- Thus we consider 
Hamiltonians of the form 

H =-j £ f i+cos2( ^-^y. ( 2 ) 



<i,j>& 
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In JTU] we call a bond "ordered" if the angle 9 between two neighboring sites 
is small enough. Here we call it ordered if the angle (9mod7r is small enough. 
Then the argument of [TOj , based on (2H] and [T7j goes through without any 
changes. There is a first-order phase transition for p chosen large enough 
(in general the values of p for which the proof works depend on N and d) 
between a high-temperature regime in which most nearest neighbor bonds 
are disordered, and a low-temperature regime in which most nearest neighbor 
bonds are ordered. This holds for each dimension at least 2, and whereas the 
Mermin- Wagner theorem excludes nematic long-range order in d = 2 in 
d = 3 and higher long-range order will occur PJ. Between the ordered and 
the disordered phase(s) free energy contours occur, whose probabilities are 
estimated to be uniformly small via a contour estimate valid over a whole 
temperature interval. In the contour estimate use is made of the Reflection 
Positivity [22j of the model. We remark that at low temperatures percolation 
of ordered bonds holds jHj; it follows from our results that the associated 
percolation transition is first order. 

For lattice gauge models the variables are elements of a unitary repre- 
sentation of a compact continuous gauge group, e.g. U(l), SU(n), or sums 
thereof [30 . Here we present the argument in the simplest case of a U(l)- 
invariant interaction in 3 dimensions. 

H = -J ( 3 ) 

plaquettes PsZ 3 
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with L(U P ) = ( 1+cos (^i+^ ^ 3 ^) Y Here the e, denote the 4 edges 



making up the plaquette P. 

The effect of choosing the nonlinearity parameter p high is that the po- 
tential, although it still has quadratic minima, becomes much steeper and 
narrower. In this way one constructs in a certain sense a "free energy barrier" 
between ordered and disordered phases. 

The lattice gauge model proof becomes similar to the arguments from 
|17j . When the product over the link variables is sufficiently close to unity, 
we'll call the plaquette "ordered", "disordered" otherwise. This distinction 
corresponds to unfrustrated and frustrated plaquettes in the Potts case. We 
will sketch the argument for the toy model where the potential L(U) is chosen 
to be 1 if the sum <\>p of the oriented angles along the plaquette P is between 
— | and +| and zero otherwise. The generalization to the non-square well 
potentials can then be done in the same way as in ^U|. The correspondence, 
as in jTUj, is that e is of order 0(-^=). 

Our strategy is to find bounds for free energy contours between ordered 
phases, in which one has mainly cubes with 6 ordered plaquettes, and dis- 
ordered phases, in which most cubes have 6 disordered plaquettes. We need 
thus to estimate the weights of contours consisting of cubes which are neither 
ordered nor disordered. The number of possibilities for such cubes includes 
the 7 possibilities given in [T7| . except that now we have the additional 8th 
possibility of having cubes with one disordered plaquette on one side, and 
five ordered ones. 

For the partition function Z L on a cube B L of size L 3 we use the (quite 
rough) lower bound 

Z L >max(ZiZ° L ), (4) 

where Zf (resp., Z£) is part of the partition function, calculated over all con- 
figurations which have all plaquettes disordered (resp., mostly ordered). For 
the disordered partition function Zf we obtain the lower bound (1 — 4e) 3L 
(we take a normalized reference measure, giving a weight 1 to each link). 

For the ordered partition function Z° L we proceed as follows: we first 
choose a set of bonds Tl in B>l, which is a tree, passing through every site. 
For example, we can put into Tl all vertical bonds - z-bonds - except these 
connecting sites with ^-coordinates and 1, plus all ?/-bonds in the plane 
z — 0, except these connecting the sites with y-coordinates and 1, plus 
all x-bonds of the line y = z = 0, except the one between the sites (0,0,0) 
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and (1,0,0) . The site (0,0,0) can be taken as a root of T L . Note that the 
number of bonds in Tl is L 3 — 1. Therefore it is not surprising (and easy to 
see) that for every edge configuration <fi = {(fib, b G Tl} there exists a unique 
site configuration ip = ifi^ = {ifi x ,x G Bl} , such that the following holds: 

1. Let 0^ denote the gauge transformation, defined by the configuration 
ip. Then [q^ o 0) =1 for every bond b G Tl; 

b 

2- ^(0,0,0) = 1- 

For every family of bonds S C B L let us define a bigger family C (S) , by the 
rules: 

1. ScC(S), 

2. for every four bonds {b±, ...,64} , making a plaquette, such that three 
of them are in S, we have {bi, 64} C C (S) . 

Then we can consider also the sets C 2 (S) — C (C (S*)) , C 3 (S*) , and so on. 
Define C(S) = U^C* (S) . Note that the number of plaquettes in (£(Xl) is 
3L 3 — O (L 2 ) . We claim now that for every configuration <fi>T L = {(fib, b G Tl} 
one can specify (in a continuous way) a collection of arcs 
{I b = h (0tJ C S 1 , 6 G £ (T L ) \ T L , \I b \ = |} , such that for every configu- 
ration (fi on Bl, which coincides with (fix L on Tl, and for which the values 
(fib on the bonds b G \ Tl belong to the above segments lb, all the 

plaquettes that fall into £ (Tl) are non-frustrated. That would imply that 

2L 3 



Z° L >Q exp{3J(L 3 -0(L 2 ))} 



by Fubini's theorem. To see the validity of our claim, consider first the 
case when the configuration <fi>T L = 1 G S 1 (here 1 is the neutral element). 
Then the choice of the segments I b is easy: i&(l) = [— §, §] for every 
b G €(T L ) \ T L . For a general (fi TL let us take the corresponding gauge 
transformation q^ t l (which is the identity for 4>t l = 1), and we define our 
segments by 

This provides a lower bound 



Zl > max 



l-Asf L3 ,Q 2L exp{3J(L 3 -0(L 2 ))} 



(5) 
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This bound on the partition function as the maximum of the ordered and 
disordered term is similar to the argument in ^U|. It plays the same role as 
the bound in terms of a fixed energy partition function given in [T7] . 

To obtain our contour estimates, by Reflection Positivity we need to esti- 
mate the partition functions of configurations constrained to have a "univer- 
sal contour". The estimates of the 7 types of universal contours mentioned 
in [T7] are of a similar form as in that paper with the number of Potts states 
q up to some constants replaced by -. The universal contour due to the 
new case of cubes with one disordered plaquette consists of configurations 
in which the horizontal plaquettes in every other plane are disordered, and 
all the other ones are ordered. These configurations have a similar entropy 
contribution to the partition function as the ordered configurations, but the 
energy per cube is | of that of a cube in the fully ordered situation. For 
e small enough (which corresponds to p large enough) also such a contour 
is suppressed exponentially in the volume. The combinatorial factor in the 
contour estimate changes by some finite constant, which choosing e small 
enough takes care of. 

To summarize we have obtained the following results: 

Theorem 1 For any nonlinear RP N ~ l model in dimension 2 or more and 
p high enough, there is a first order transition, that is, there exists a temper- 
ature at which the free energy is not differentiable in the temperature param- 
eter. 

Theorem 2 For lattice gauge models with plaquette action (^-^^j , (where 

L(Up) = Tr(Up + Up)) in dimension 3 and more and p high enough, there 
is a first order transition, that is there exists a temperature at which the free 
energy is not differentiable in the temperature parameter. 

Here Up denotes the adjoint operator of Up. 

This seems to be the first case in which a first order transition for a lattice 
gauge model with a continuous gauge symmetry group is rigorously obtained. 
Whereas the example of the Potts lattice gauge model in d = 3 or higher is 
between a confining and a nonconfining phase fHl E] m our theorem this 
is to be expected in d — 4, with U(l) symmetry only. For d = 3 and also for 
SU(n) in d = 4 we conjecture that confined phases exist on both sides of the 
phase transition. 

Our proof only gives results for very high p. We will discuss some further 
aspects of what may actually be the p- values for which to expect first-order 
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transitions, and what one might hope to prove. The recent work of Biskup 
and Chayes shows that if a reflection positive model has a phase transition 
in mean field theory, then also at sufficiently high dimension a first-order 
transition occurs. They include in their discussion the RP N ~ Y model for 
N = 3, for which even for p — 1, a first-order transition is derived. The 
mean field analysis of j3] indicates that a similar result for the ferromagnetic 
case holds if p = 3. For lattice gauge models also the standard actions lead to 
first order transitions in mean-field theory [33], section 34.4, which indicates 
a first order transition in sufficiently high dimension. 

If one believes that here the spherical (N to infinity) limit is not singular 
(which has been a matter of controversy itself), then for the square lattice, iV 
large and p larger than 6 the ferromagnet might have a first order transition, 
while for the RP 1 ^^ 1 case on the square or triangular lattice even for p = 1 a 
first order transition occurs, while for the hexagonal lattice one presumably 
needs a higher value of p jSH EH] ■ 

As mentioned before, numerical work suggests that in the standard (p = 
1) Lebwohl-Lasher model with N = 3 as well as in the [/(l)-lattice-gauge 
model a first order transition should occur; however, this seems far away 
from any provable result. 
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